Trident pair production in strong laser pulses 
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We calculate the trident pair production amplitude in a strong laser background. We allow for finite pulse dura- 
tion, while still treating the laser fields nonperturbatively in strong-field QED. Our approach reveals explicitly 
the individual contributions of the one-step and two-step processes. We also expose the role gauge invariance 
plays in the amplitudes and discuss the relation between our results and the optical theorem. 
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Electron-positron pair production, at the focus of an intense 
laser, is currently a topic of considerable interest due to the 
development of extreme light sources such as ELI (]]]. As 
is typical for particle scattering experiments, many different 
processes may contribute to the final yield of pairs. It is im- 
portant to be able to judge the relative importance of these 
processes, and distinguish their contributions from each other. 
In the SLAC El 44 experiment |01, pairs were produced by 
colliding a (low intensity) laser with the SLAC electron beam; 
high energy photons radiated by the electrons then combined 
with photons in the laser to produce pairs. In that experiment, 
a process called 'trident' may also have produced pairs. How- 
ever, its contribution could only be evaluated approximately, 
since no exact expression for the trident amplitude was avail- 
able, nor has one been given to date Jjl] . Accounting theoret- 
ically for all relevant effects in modern laser experiments is 
a formidable challenge. As well as contributions from pro- 
cesses such as trident (along with, for example, vacuum pair 
production J3l and cascades J5t]), one would like to include 
effects due to the properties of the laser itself, such as their 
ultra-high intensity (currently 10 22 W/cm 2 ), ultra-short dura- 
tion (measured in femtoseconds or even attoseconds) and tight 
focus (focal diameter of the order of 10 microns). 

In this Letter, we present the first complete calculation of 
the trident process. We include high-intensity effects by treat- 
ing the laser field nonperturbatively and exactly. We include 
finite size effects by allowing finite pulse duration. We be- 
gin by describing the essential features of the trident process, 
modelling the laser as a null field and using Volkov solutions 
to describe the scattered particles. Key aspects of existing cal- 
culations are discussed, and shown to be unphysical due to vi- 
olations of the Ward identities of QED. We then calculate the 
scattering amplitude, emphasising its physical content, and re- 
late our results to the optical theorem. We use gauge invari- 
ance to simplify the final expression for the emission rate, and 
discuss approximations useful to the high intensity regime. 

We consider an electron, incident upon a laser field, emit- 
ting a photon. This photon then combines with photons in 
the laser to produce an electron-positron pair. The relevant 
Feynman diagram is shown in Fig. Q] Double lines repre- 
sent the fermion propagator in the background laser field. We 
take this background to be a plane wave, with field strength 
Fpv = F )iv (k-x), where the momentum is lightlike, so 
k 2 = 0. The laser frequency is (0 = |k|. Our treatment holds 
for arbitrary k-x dependence: in particular, we have finite pulse 
duration when F^y vanishes, or goes rapidly to zero, outside 
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FIG. 1: The trident process. Double lines represent the fermion prop- 
agator in the background field. The second term accounts for the 
exchange of indistinguishable fermions. 



of some k-x range. 

Fig. Q] actually describes two processes. The first is 'one 
step', in which the intermediate photon is virtual: this is the 
process traditionally referred to as trident. The second pro- 
cess is 'two step', in which a real photon is scattered from the 
incoming electron (allowed in a background field via nonlin- 
ear Compton scattering ffj), and this real photon then goes 
on to create a pair via stimulated pair production |7|]. This 
was the process of interest in the SLAC experiment. Since 
the Feynman diagram in Fig. Q] makes no distinction between 
these processes, we will refer to the full diagram as 'trident', 
and use the notions of one- and two-step processes to distin- 
guish the two contributions. The trident S-Matrix element is 
easily written down, using strong-field QED, in terms of the 
photon propagator G^y and Volkov solutions \jfj carrying the 
asymptotic momenta pj in Fig.[T|ii : 
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The Volkov solutions also carry a dependence on k-x, and k-y, 
which describes the "dressing" of the particles by the back- 
ground. One trades this dependence, via Fourier transform, 
for new variables r and s: in the periodic plane wave case, 
the dressing is responsible for the shifted electron mass, and 
the Fourier transform reduces the S-Matrix element to a dis- 
crete sum over processes involving different numbers of 'laser 
photons' i^llTofl. Our continuous r and s are the analogues of 
photon number for pulsed plane waves, see 0]. After this, 
the x and y integrals trivially yield two delta functions con- 
serving four-momentum. One of these fixes the intermediate 
photon momentum k! , and we have, writing 8p = p\ — pi and 
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p out = P2+ Pi + Pa from here on, 

S fi = (27t) 2 e 2 jdrjds M(r,s)8 4 (p oul - p x -rk- sk) 
1 
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The expression for M will be given soon: it is lengthy and 
we do not need it yet. The Fourier transformed (Feynman 
gauge) photon propagator is g^v/ik' 2 + is). The remaining 
delta function, above, sees the asymptotic momenta and hence 
the free electron mass, not the shifted mass of the infinite 
plane wave calculations; this is a consequence of allowing fi- 
nite pulse duration, just as is observed in 10] for stimulated 
pair production and in Ml ill for nonlinear Compton scattering 
in pulsed fields. From the delta function in (0 it becomes 
clear that r and s parameterise the energy-momentum taken 
from the laser. 

In the current literature, a divergence has been identified 
in (0, arising from the pole in the photon propagator in the 
two-step process, i.e. when enough energy is taken from the 
background to put the photon on-shell and k' 2 = 0. This diver- 
gence was attributed to the infinite temporo-spatial extent of 
the (periodic plane wave) background considered. Based on 
this, it was suggested that the divergence could be dealt with 
by modifying the photon propagator according to IU2I1 
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which damps the propagator outside of a chosen, fixed, tempo- 
ral range T, in order to model finite pulse duration. However, 
it is well known that the structure of propagators and vertices 
in QED is severely restricted by gauge invariance, violations 
of which are measured by the Ward identities fl3ll . The pre- 
scription (0 manifestly violates the Ward identity 11411 



(4) 



which is sufficient to reject (0. We briefly state some addi- 
tional physical reasons. First, the origin of the divergence can- 
not be the infinite extent of the background. Any background 
which permits the two-step process, producing an on-shell in- 
termediate photon, will admit the k' 2 = divergence, indepen- 
dent of the background's spacetime support (even if the theory 
contains a natural analogue of the cutoff scale T, for example). 
Secondly, there is also no reason why any chosen T should 
be preferred by the theory over any other. The prescription 
(0 also breaks manifest Lorentz invariance. While this may 
not be considered too serious, given that the background field 
already introduces preferred directions in spacetime, these di- 
rections are lightlike, not timelike as 01 implies. Using 01 
also gives a T dependent correction to the divergence-free 
one-step process, and leaves unphysical factors of T in the 
weak field limit where all background effects should become 
negligible. It is therefore difficult to ascribe a physical mean- 
ing to results following from 01. 



The resolution of the above problem is simple: in a com- 
plete treatment, there is no divergence. It appears only when 
one neglects the pole prescription already contained in Feyn- 
man propagators. The photon propagator is really 1 / [k' 2 + is) 
and one must remember to take the infinitesimal e — > + . To 
show there is no divergence, we return to 01, where the prop- 
agator is evaluated at k! = dp + rk. This can certainly go on- 
shell, so k' 2 = 0, for certain r. This behaviour may be exposed 
by recalling the distributional result: 
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which clearly separates contributions from, respectively, real 
(on-shell) and virtual (off-shell) photons. We will reaffirm 
this statement below. Here, we perform the essential step, 
inserting (0 into 0). In the first term, the delta-function is 
8(k' 2 ) = 8(2rk-8p + 8p 2 ) and is therefore eliminated by per- 
forming the r integral. This fixes 
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so the propagator's pole fixes, quite naturally, the energy 
transfer from the laser such that the photon goes on-shell, 
i.e. such that the two-step process occurs. There is no di- 
vergence here, nor any need to deform the theory along the 
lines of 0). (Note, the denominator of r 1 is nonzero because 
k-8p = k- (pi + Pa) > 0, using the momentum conservation 
law in (0.) To proceed, one eliminates the s integral using 
one component of the delta function in (0. This is easily 
done in lightfront co-ordinates x ± := x° ±x 3 , x L := {x 1 ,* 2 }, 
and q± = (qo±qT,)/2 for momenta. Since the laser momen- 
tum kfj is lightlike, we may choose k + to be its only nonzero 
component. Hence, s appears only in the q + component of the 
delta function, which eliminates the s integral and sets 
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Three components of the delta function remain, conserving 
momentum in the — and _L directions, transverse to the laser. 
Writing 8"(q) = S 2 (q ± )S(q-)/k + , the S-Matrix element is 
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This completes the calculation of 5//. It should be clear from 
(0 that the off/on-shell parts of (0 correspond precisely to 
the one/two-step processes. Nevertheless, let us confirm this, 
which will also serve as a check on Sfi before moving on to 
the cross section. Consider first the one-step process, with 
overall momentum conservation given by the explicit delta- 
function in (0. Squaring the argument of the delta function, 
one finds that r + s > 4m 2 /(k-pi), which states that the to- 
tal incoming energy (that of the initial electron and that taken 
from the laser), must be sufficient to produce three particles 
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of rest mass m. This is the only constraint in 5/,- if the pho- 
ton is off-shell, k' 2 ^ 0, and we pick up the principal value 
term. This is the 'trident' contribution in the old nomencla- 
ture. (Again, there is no dependence on the shifted mass found 
in a periodic plane wave.) Consider now the two-step process, 
which has two momentum conservation relations 



P\ + rk = p2 + k' nonlinear Compton scattering 
k' + sk = pj,+ P4 stimulated pair production . 
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Squaring, one finds that, for the individual processes to occur, 
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Using the overall delta function in (0, it is straightforward to 
show that the fixed parameters r 1 and s,j which appear in our 
on-shell part obey r 3 > and r s i > 2m 2 /k-8p. Recalling that 
k-8p = k-k' evaluated on shell, we recover dTOb . Our solution 
thus contains the correct kinematics of both the one- and two- 
step processes, which are described precisely by the off- and 
on-shell parts of the decomposition (0. Thus, our solution 
allows the contributions of these processes to be individually 
calculated and compared. 

The identity (0 also corresponds to a split into real and 
imaginary parts. From the optical theorem for scattering am- 
plitudes (that is, unitarity of the S-Matrix), one expects the 
appearance of imaginary, or absorptive, parts to correspond 
to the excitation of real rather than virtual intermediate states 
[16, [l7ll . Indeed, we have seen that the imaginary part of 
corresponds to the intermediate photon becoming real. 
The existence of this imaginary part is entirely due to the 
dressing of the fermions by the background: in cutting lan- 
guage, the dressing allows one to perform a cut through Fig.Q] 
which yields physical, non-zero scattering amplitudes; those 
for nonlinear Compton scattering and stimulated pair produc- 
tion. (The trident amplitude itself may be obtained by cutting 
the two-loop fermion propagator IU8I1 .) Our result therefore 
appears to contain a rather novel example of the optical theo- 
rem at tree level, made possible by the background. We hope 
to investigate this further in the future. 

In the remainder of this Letter we give the complete cal- 
culation of the emission rate. We will also show how gauge 
invariance simplifies the final results. 

Our background field has potential A )1 (k-x) = fj(k-x)a^, 
where j is summed over the transverse directions. The po- 
larisation vectors obey a' -k = 0, a' = —m 2 a 2 /e 2 8^ which 
defines an invariant, dimensionless amplitude a. From the cor- 
responding Volkov solutions, we define 
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These are combined into 'nonlinear Compton', F^, and 'pair 
production', A^, parts as follows (we denote the reverse- 



ordered S by S): 

r^(k-x):=u P2 S{p2,k-x)fS{p h k-x)u Pi 
x exp{U(p2,k-x 7 °°) + U(pi, 
^(k-y):=u P3 §{p 3 ,k-y)^S(-p4,k-y) Vl 
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The limits in J are prescribed by the LSZ reduction formula 
for the trident amplitude. The product of the Fourier trans- 
forms of these parts gives us the amplitude M(r, s) from above, 
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To expose the further role of gauge invariance in our result, 
we give the explicit form of the pair production part A' 1 (from 
which one may write down by crossing symmetry): 
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All dependence on the parameter s (a 'photon number' ) is con- 
tained in four functions B0...B3. For j = 1,2,3, and writing 
/3 = /j + f 2 , these are defined by (sum over n = 1,2,3) 



Bj(s):= /dz/ ; (z)exp 



isz + ia„ / dw /„ (w) 
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where the coefficients a„ may be read off from (fTTT) -([T3l): 
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The functions dT5l > are finite, provided the fj vanish asymp- 
totically, as is the case for pulsed fields. We assume this be- 
haviour from here on. (The periodic plane wave case may be 
recovered in a suitable limit, see lfl5ll .) The fourth function Bq 
is defined as in ( TT31 ) but without any damping factor of / un- 
der the integral. Consequently, Bq is the Fourier transform of 
a pure phase, and some prescription is required for calculat- 
ing it, at least numerically. Regularisations of this integral (or 
its equivalent in the nonlinear Compton part of the amplitude) 



were suggested in fll9l 12011 . but there is a more fundamental 
method of defining Bq : one appeals to gauge invariance. Mak- 
ing a (quantum) gauge transformation in (0, one finds that Sfj 
is gauge invariant provided that (sum over n = 1,2,3,) 

sB (s) = O n B n (s), (17) 

and similarly for analogous functions in . Hence, by defin- 
ing Bq in terms of the well-behaved By, gauge invariance re- 
duces the number of functions in play from four to three. This 
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result should be compared with equations (A3) in |@], (6) in 
fioll . and (23) in 1 11 1: those expressions were obtained either 
as identities of special functions or by regularisation, but they 
are unified and explained physically by gauge invariance. 

One now has all the necessary ingredients to calculate ((8]). 
To obtain the full rate, one takes |S/;| 2 and divides out the vol- 
ume transverse to the pulse, as usual. The remaining steps 
are to sum/average over spins and perform the final state inte- 
grations, where again it is natural to use lightfront variables. 
There are nine integrals, corresponding to three momentum 
components for three outgoing particles. Three of these are 
eliminated by the remaining 5 lf . Writing K for the bracketed 
term in ((8), one obtains, with a = e~ /(All), 
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The instruction 'shell' indicates that each p + is evaluated on 
shell, i.e. p + = (p\ + m 2 )/(Ap_), and p2, which is eliminated 
by momentum conservation, obeys p2 X = (p\ — P3 — P4)i 
and p2_ = (pi — P3 —Pa)-- The remaining integrals are over 
the momenta of the produced electron-positron pair. The full 
emission probability is a function of the laser amplitude, a, 
the pulse geometry, and the incoming momenta which appear 
through k-p\. This completes the calculation. 

Despite the length of the expressions involved, our final re- 
sult for the rate is not more complicated than that of the pe- 
riodic plane wave case studied in 11211 . Let us compare the 
expressions therein with our own results. There are three real 
differences. The first is that the discrete sums over photon 
number are replaced by Fourier integrals over r and s. The 
second is that the shifted mass plays, in general, no role in our 
expressions. Finally, our S-Matrix contains two distinct terms 
corresponding to the one- and two-step processes. 

An important step in the calculation of the rate is the evalua- 
tion of the functions dT3T >. In general, these must be calculated 
numerically. In various limits, though, analytical approxima- 
tions exist. The weak field limit is an obvious example: one 
simply expands everything in powers of the field amplitude a. 
What is more interesting for modern experiments is the high 
intensity limit, which amounts to a 3> 1. This is studied in 
I IM m the context of nonlinear Compton scattering, using an 
asymptotic expansion of the (equivalents of) the Bj. Let us 
adapt this to our trident process: the idea is to look for points 
of stationary phase in (fT~5T >. However, one can show from the 
kinematics that s — ajfj(k-x) ^ QVk-x and hence no points 
of stationary phase exist. Thus, one can try to instead deform 
the contours into the complex plane and use the method of 
steepest descent. This is a well understood technique and may 
be applied immediately to the trident amplitude. It would be 
interesting to try and understand the physics behind this ap- 
proximation in more detail: as a first step, one could return to 
the infinite plane wave case and attempt to establish a connec- 
tion between the saddle points in the complex plane and the 
photon number, or the effective electron mass. 



In conclusion, we have given the first full calculation of tri- 
dent pair production in a laser field, using strong-field QED. 
We have included finite size effects due to the ultra-short du- 
ration of modern pulses. Both the one-step and two-step pro- 
cesses involved have been explicitly identified, and our results 
are in agreement with the optical theorem. We have also re- 
vealed the role gauge invariance plays, not only in the trident 
process, but also in nonlinear Compton scattering and stimu- 
lated pair production, through the relation ( TPTl i. We remark 
that our approach is equally valid for laser assisted M0ller 
scattering, as this is just the crossed process of trident pair pro- 
duction: the appropriate S-Matrix element is obtained from 
(0 by taking the outgoing positron to be an incoming elec- 
tron. 

Contrary to recent claims in the literature, we have shown 
that there is no divergence in the trident amplitude. Despite 
this, the numerical methods previously employed to calculate 
the amplitude, from which a great deal of information was 
obtained, are equally appropriate here. In particular, one may 
now easily compare the contributions of the one-step and two- 
step processes, as well as test the old Weizsacker- Williams ap- 
proximation for the former J3| ■ It will be extremely interesting 
to see what these investigations reveal. 
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